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Abstract 

We assume the t-t'-J model to describe the Cu02 planes of hole-doped cuprates 
and we adapt the spin-charge gauge approach, previously developed for the t-J 
model, to describe the holes in terms of a spinless fermion carrying the charge 
(holon) and a neutral boson carrying spin 1/2 (spinon), coupled by a slave-particle 
gauge field. In this framework we consider the effects of a finite density of incoherent 
holon pairs in the normal state. Below a crossover temperature, identified as the 
experimental "upper pseudogap", the scattering of the "quanta" of the phase of the 
holon-pair field against holons reproduces the phenomenology of nodal Fermi arcs 
coexisting with gap in the antinodal region. We thus obtain a microscopic derivation 
of the main features of the hole spectra due to pseudogap. This result is obtained 
through a holon Green function which follows naturally from the formalism and 
analytically interpolates between a Fermi liquid-like and a d-wave superconductor 
behaviour as the coherence length of the holon pair order parameter increases. 
By inserting the gauge coupling with the spinon we construct explicitly the hole 
Green function and calculate its spectral weight and the corresponding density of 
states. So we prove that the formation of holon pairs induces a depletion of states 
on the hole Fermi surface. We compare our results with ARPES and tunneling 
experimental data. In our approach the hole preserves a finite Fermi surface until the 
superconducting transition, where it reduces to four nodes. Therefore we propose 
that the gap seen in the normal phase of cuprates is due to the thermal broadening of 
the SC-like peaks masking the Fermi-liquid peak in the spectral weight. The Fermi 
arcs then correspond to the region of the Fermi surface where the Fermi-liquid peak 
is unmasked. 



PACS numbers: 71.10.Hf, 11.15.-q, 74.72.-h, 74.72.Kf 



1 Introduction 



The phenomenon of "pseudogap" in hole-doped cuprates appears rather complex, exhibit- 
ing different features in various regions of the doping-temperature phase diagram and even 
the word "pseudogap" is used with different meanings by different authors. One of the 
most spectacular features associated with pseudogap phenomenology is the appearance 
revealed by ARPES of a truncated Fermi surface consisting of Fermi arcs. There have 
been several theoretical explanations of this phenomenon, many of them critically ana- 
lyzed in [1] , but no consensus has been reached, since none of these proposals explains all 
the characteristic features of the experimental data. 

In this paper we develop an "explanation" of the Fermi arcs within a generalization 
of a gauge approach to superconductivity in cuprates recently proposed [2], comparing 
the results with ARPES and tunneling data. To understand our proposal it is useful 
first to sketch the pairing mechanism, eventually leading to superconductivity, for under- 
doped cuprates presented in [2] within a gauge approach to the t-J model: as we dope a 
vortex-like quantum distortion of the AF background is generated around the empty sites 
(described in terms of fermionic spinless holons) with opposite chirality for cores on the 
two Neel sublattices. The spin excitations (bosonic spin- 1/2 spinous) are gapless without 
doping, corresponding to long-range AF order, but above a critical doping density they 
acquire a finite gap due to scattering against the spin vortices and the long-range anti- 
ferromagnetic order is converted to a short-range order. Due to the no-double occupation 
constraint, decomposing the hole into holon and spinon generates a local gauge symme- 
try inducing in turn a gauge attraction between holon and spinon binding them into a 
physical hole. 

The primary pairing force for the charge carriers is an attraction due to chirality be- 
tween spin vortices with cores on two different Neel sublattices, inducing an attraction 
between holons in their cores. As a consequence of this attraction at a crossover tem- 
perature, denoted Tph, a finite density of incoherent holon pairs are formed. In [2] it 
was claimed that this phenomenon produces a reduction of the hole spectral weight on 
the Fermi surface (FS) and it was proposed to identify this temperature with the exper- 
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imentally observed "upper pseudogap" , where the in-plane resistivity deviates from the 
hnear behaviour. At a lower crossover temperature, denoted Tpg, also a finite density of 
incoherent spinon pairs appears, giving rise to a gas of incoherent preformed hole pairs 
through holon-spinon gauge attraction . 

Finally, at a even lower temperature, Tc, the hole pairs become coherent and super- 
conductivity emerges. This approach exhibits yet another crossover temperature [3], T*, 
intersecting Tpg in the doping-temperature phase diagram. Such crossover is not directly 
related to superconductivity. It corresponds to a change in the holon dispersion. It 
is characterized by the emergence of a "small" holon Fermi surface around the momenta 
(±7r/2, ±7r/2), with complete suppression of the (coherent) spectral weight for holes in the 
antinodal region and partial suppression outside of the magnetic Brillouin zone (MBZ). 

This crossover appears only in bipartite lattices. Below T* the effect of short-range 
AF fluctuations becomes stronger and the transport physics of the corresponding normal 
state region is dominated by the interplay between the short-range anti-ferromagnetism 
of spinous and the thermal diffusion induced by the gauge fluctuations triggered by the 
Reizer [4] momentum. This interplay produces in turn the metal- insulator crossover [3]. 
We identify T* in experimental data with the inflection point of in-plane resistivity and 
the broad peak in the speciflc heat coefficient [5]. The region above T* in the doping(5)- 
temperature(T) phase diagram will be called the "strange metal phase" (SM), the one 
below will be called the "pseudogap phase" (PG). Actually there is no agreement between 
the experimentalists on the existence of two crossovers (in our approach Tph and T*) asso- 
ciated to the pseudogap phenomenology and the same is true for most of the theoretical 
approaches. For clarity in Fig. [1] we present a schematic behaviour of the crossovers 
T*,Tph and Tpg in the phase diagram. 

In this paper we extend the pairing mechanism for holes developed for the PG of the 
t-J model to the normal phase of the t-t'-J model showing explicitely that T* and Tph are 
indeed distinct. 

It turns out that a finite density of incoherent holon pairs induces an angle dependent 
reduction of the (physical) hole spectral weight on the FS, starting from the antinodal 
regions, thus producing a microscopic derivation of the phenomenology of Fermi arcs 
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coexisting with gap in the antinodal region. Notice that this effect on the holes occurs 
even if the pairing is only among unphysical holons, the spinons being still unpaired, so 
that there is not yet a gas of preformed hole-pairs. 

We show that, when a gas of incoherent ho Ion pairs is present, an energy scale sepa- 
rating low energy modes with a Fermi liquid (FL) behaviour from high energy modes with 
a d-wave superconducting behaviour results naturally and self consistently. This energy 
scale will be identified with the the inverse correlation length ("mass") of the quanta 
of the phase of the holon pairs field. The scattering of these excitations against holons 
produces in the holon Green function lowering T a gradual reduction of the spectral weight 
on the FS at small frequency as we move away from the diagonals of the Brillouin zone. 
Simultaneously at larger frequencies we have the formation and increase of two peaks of 
intensity precursors of the excitations in the superconducting (SC) phase. The value of 
in fact decreases with T , in an extended range approximately linearly, and this decrease 
drives the system towards the superconducting phase occurring when, possibly disconti- 
nously (as experimentally suggested e.g. by [6]), jumps to 0, in correspondence with 
the condensation of the hole pairs. Decreasing m^, the well defined FL quasi-particles 
appearing at high T gradually lose their coherence in favor of SC-hke excitations. The 
latter gain spectral weight and become increasingly well defined excitations. In the SC 
phase, all the modes are above and the holon system is a (i-wave superconductor, in 
particular the holon spectral weight at the FS is reduced to zero except on four nodes. 

The physical hole is obtained as a holon-spinon resonance produced by the gauge 
attraction and it inherits the above holon features, but with a strongly enhanced scattering 
rate, due to the spinon contribution. 

The behaviour of the spectral weight derived from the above sketched mechanism 
appears consistent with many ARPES data and it is able to explain the evolution of 
the density of states as derived in tunneling experiments. This mechanism of spectral 
weight suppression exhibits the fingerprint of the presence of the slave-particle gauge 
field, because the smooth interpolation between FL and SC behaviour is actually due 
to the interaction of the phase of the holon pairs with the gauge field. Without such 
interaction the SC-like peaks are strongly suppressed outside of the SC phase, disagreeing 
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with experiments. 

Although many phenomenological features in our approach are similar to those of the 
approach proposed phenomenologically in [7] and partially justified microscopically in [8], 
at odds with those approaches in ours the hole for 7^ (i.e except in the SC phase) has 
always a FS without gap. This is consistent with the fact that the "pseudogap" region 
can be reached (at least in our approach) from a FL behaviour without crossing a phase 
transition, only crossovers being involved in the process, which appears to agree with 
experiments. We propose therefore that the gap seen in the normal phase of cuprates is 
due to the thermal broadening of the SC-like peaks, the Fermi arcs corresponding to the 
region of the Fermi surface where the Fermi-liquid peak is unmasked from that broadening. 
A further complete suppression of the spectral weight in the antinodal region occurs in 
the PG, as previously discussed. Readers only interested in comparison of final theoretical 
results with experiments can skip the next sections and go directly to section 8, where we 
recapitulate the key points of our approach before comparing with experimental data. The 
effects of the reduction of the spectral weight on transport properties will be discussed in 
a companion paper [9] where we prove, in particular, a deviation with negative curvature 
from linearity in T of in-plane resistivity. 

2 The spin-charge gauge approach to the t-t'-J model 

The approach in its original version [3] assumed the t-J model as model Hamiltonian 
for the low-energy physics of the CUO2 planes of the cuprates. Here we add a negative 
ncxt-ncarest-neighbor hopping term t' (with \t'\ < t) which is known necessary to fit the 
FS of cuprates. The scheme of the approach, however, remains unchanged and we sketch 
it here for reader's convenience and to set up the notation, emphasizing only the new 
features. 

We decompose the hole operator Cjo- at site i of the t — t' — J model as Cia = h*bia, 
where h is a spinless fermionic holon, carrying charge, while is a spin 1/2 bosonic 
spinon, carrying spin and obeying the constraint ^^b*^bia — 1. The redundant degree 
of freedom arising from this decomposition is cured by an emergent slave-particle U{1) 
gauge field, A^, minimally coupled to holon and spinon with the same charge. With the 
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choice of statistics adopted the no-double occupation constraint is automatically satisfied 
because of the Pauli principle for the holon. [ For simplicity in this paper we use the same 
symbols to denote the field operators in the hamiltonian formalism and the corresponding 
fields in the path-integral lagrangian formalism. ] One then uses the possibility offered 
in 2D (and ID) to add a "statistical" spin fiux (e**") to ba and a "statistical" charge 
flux (e"'*'') to h "compensating" each other so that the product e''^^^he^^% is still a 
fermion. The introduction of these fluxes in the lagrangian formalism is materialized 
through Chern-Simons gauge fields. We then optimize their choice in an improved mean- 
field approximation (MFA). A key step of MFA is to find a reference spinon configuration 
with respect to which expand the fiuctuations that will be described by a new staggered 
spinon field, Za- This reference configuration is found optimizing the free energy of holons 
in the presence of a fixed holon-dependent spinon configuration. One can show as in 
the ID [10] and 2D [11] t-J model that the optimization involves a spin-fiip associated 
to every holon jump between different Neel sublattices. Furthermore, if we neglect the 
spinon fluctuation Za in the spin flux, the effect of the optimal spin flux is to attach a spin- 
vortex to the holon, with opposite chirality on the two Neel sublattices. More precisely 
with this approximation 



The gradient of can be seen as the potential of a vortex. These vortices take into 
account the long-range quantum distortion of the AF background caused by the insertion 
of a dopant hole, as flrst discussed in [12], the rigidity holding up them being provided 
precisely by the AF background. 

Neglecting the holon fluctuations in the optimal charge flux in the 2D t-J model 
was argued to provide a staggered flux tt per plaquette for small enough doping and 
temperature. This anomalous behaviour near half-filling was justified on the basis of a 
rigorous result by Lieb [13] and a numerical simulation, proving that in a square 2D lattice 
at half-filling the optimizing flux for the free energy of fermions is indeed translationally 
invariant and tt per plaquette. The contribution from the reference spinon conflguration 
in the optimization can be shown to trivialize the optimal flux at high enough doping 




(2.1) 
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or temperature, in agreement with diamagnetic inequality. Although Lieb's proof does 
not extend to the t-t' model, one can still rely on the results of [14] which show that 
for elementary circuits, triangles and squares, the optimal flux at half- filling is ±7r/2 for 
triangles and ivr for squares. Assuming that the above rule for fluxes matches with the 
staggered flux vr per square plaquette of the underlying t model one founds that along the 
diagonal links of the t-t' model the flux should be trivial, this implying that parity and 
time-reversal are still preserved. On the basis of the above considerations we extend the 
conjecture made for the t-J model as 

Assumption The optimal charge flux in the t-t'- J model at sufficiently small doping and 
temperature is given as follows: let i be a site in the even Neel sublattice and i±a, a = 1,2 
be its n.n. sites in the 1 and 2 directions respectively and i ± /3, /3 = ±i ± 2 its n.n.n. 
sites. Then 

$''(«±d) - <l>^(z) = (-l)"(±z7r/4), <l>^(z±/3)-<l>^(z) = 0. (2.2) 

Also for the t-t'-J model the contribution from the reference spinon conflguration in the 
optimization can be shown to trivialize the optimal flux at high enough doping or tem- 
perature. Extending a conjecture made for the t-J model we propose that the crossover 
between the two behaviours discussed above corresponds to the crossover between the 
"pseudogap phase" (PG) and the "strange metal phase" (SM) in the cuprates. Neglect- 
ing the fluctuations of the gauge fleld which can be reinserted by Peierls minimal 
substitution, the leading terms of the Hamiltonian can then be written as: 

H ^ En.n.<^j>i-i)^M,,h*h,e^^'''^-''^^ + J(l - h*K - h*h,){l - |AM,,n 

+Jhlh*h,K\RVB\% + Y.n.n.n.«^,»t'AM.-^^Khr (2.3) 

where AMij = Tr(z*e*^*^~*^^2;j)^*^(2), with (i) = 1 (resp.= *) if i is in the even (resp. odd) 
Neel sublattice, and it is a kind of Affleck-Marston spinon parameter [15] and RVBij = 
Yl ^ajsZiaZjjs is an RVB spinon singlet order parameter. The AM/RVB dichotomy in (12. 3p 
is due to the spin-flip in the optimization procedure described above . Then we use the 
following MFA for ho Ions and spinous: in the flrst term in (12. Sp we take < AMij >^ 1, 
while in the second term we replace the hole density by its average and in the normal state 
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we neglect the third term because of being higher order in doping (5). A long- wavelength 
treatment of the second term in (12. 3p leads to a spinon (CP^) non-linear sigma model 
with an additional term coming from the spin flux, 

J{l-26){V<^'fz*z, (2.4) 

where 9^<l''^(x) = £^,^9,^ ^^•(— l)l-'lA^^(x — j)h*hj with A the 2D lattice laplacian. In a 
quenched treatment of spin vortices one finds for the average < (V$^)^ >= ~ 6\ log6\, 
producing a mass gap for the spinon, consistent with AF correlation length at small 6 
derived from the neutron experiments [16]. 

In the parameter region to be compared with the PG of the cuprates the charge 
flux TT per plaquette converts the spinless holons h into Dirac fermions with small Fermi 
surfaces centered at (±7r/2, ±7r/2): ep ~ tkp , where kp ^ S is the (average) holon Fermi 
momentum and ep the Fermi energy. The holon dispersion is defined in the Magnetic 
Brillouin Zone. In the parameter region to be compared with the SM of the cuprates, 
where the optimal charge flux per plaquette is 0, we recover a standard tight-binding 
"large" FS ( ei? ~ t(l + 6)) for the holons, centered at (vr, vr). Let us briefly comment how 
these FS arise. In the t — t' — J model the density of the Gutzwiller projected holes is 
proportional to the doping 5, the holes corresponding in coordinate space to the empty 
sites. Substituting the representation of the hole fields in terms of holon and spinon 
in the lagrangian in the presence of the reference spinon configuration identified by the 
optimization sketched above, for the holons dressed by the charge flux in the SM one finds 
the 0-energy level at the position of the FS of the tight-binding unprojected model. This 
result is compatible with the holon density 5 if 2 holons can occupy the same momenta and 
the "vacuum" of the model is set at half- filled holon band, namely the holons relevant 
physically for the projected holes are those corresponding to the deviation from half- 
filling. Noticing [3] that the (fermionic) holon fields dressed by the charge flux describe 
semions (with statistics intermediate between fermions and bosons), the first property is 
guaranteed if we apply to them the Haldane-Wu statistics [17] for semions (as in the ID 
t — J model [10]) allowing precisely double-occupation in momentum space for spinless 
semions. We argue that this MF treatment of holons, although not strictly equivalent 
to the Gutzwiller projection for holes, nevertheless it is still reasonable for small holon 
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energies, close to the FS . In the same approximation in the PG one finds for the holons 
dressed by the charge flux the 0-energy level at the position of the small FS quoted above, 
the "vacuum" for holons then corresponds to the filling of the lower branch of the Dirac 
double cones with vertices at (±7r/2, ±7r/2) generated by the charge vr-fiux background, 
yielding a result characteristic of a 2D doped Mott insulator. 

Holons and spinous are coupled by the emergent gauge field A yielding overdamped 
resonances with strongly T-dependent life-time. This dependence originates from the 
fiuctuations of the transverse mode of the gauge field, dominated by the contribution of the 
gap less holons. Their Fermi surface produces an anomalous skin effect, with momentum 
scale 

Q ^ {Tklf'\ (2.5) 

the Reizer momentum [4]. One can take into account approximately an external frequency 
u » T replacing T by w in the life-time and in the Reizer momentum. 

Let us now turn to holon pairing. Spin vortices centered on holons on the two Neel 
sublattices have opposite vorticity and this produces an attraction, previously neglected 
in the MFA; this is the origin of the attractive force between holons. Physically it is due 
to the quantum distortion of the AF background caused by the holes. We include this 
effect in MFA by introducing also the term coming from the average of z*z in (12. 4p . We 
get the contribution: 

J(l - 25) < z*z > ^(-1)I'I+I-'IA-^(z - j)h*hih*hj, (2.6) 

where A is the 2D lattice laplacian. In the static approximation for holons (12. 6 p describes 
a 2D Coulomb gas with coupling constant J = J(l — 26) < z*z >, with < z*z >~ 
J cPq{(f + m^)^^ ~ (A^ + mD^^^ — nis, where A ^ 1 is a UV cutoff, and charges ±1 
depending on the Neel sublattice. For 2D Coulomb gases with the above parameters a 
pairing develops for a temperature Tph ~ J /2tx, which turns out to appear in the SM. As it 
occurs for non-weakly coupled attractive Fermi systems, below Tp^ there is a temperature 
at which the pairs condense, and it will turn out to be the superconducting transition 
temperature Tc, see [2]; in-between there is a crossover, Tp^, due to spinous, but it will 
not be discussed here. This omission affects only the lower temperature region for the 
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hole Green function discussed in this paper; in principle the approach is able to derive 
this correction, but it would highly complicate the calculations without affecting most of 
the main features. We defer this issue to a future publication. 

3 Low energy Hamiltonian for holons 

In this section we discuss in detail the pairing among holons in the SM region of the 
t — t' — J model in the BCS approximation, adapting the framework developed in [2], in 
turn inspired by [18]. A brief comment on the modifications needed for the PG region are 
added at the end. In the next section we incorporate in this framework the fluctuations 
of the phase of the holon-pairs order parameter. 
We start treating the kinetic hamiltonian 

HI^ = -t ^{h*hj + h.c.) -t'Y^ {h*hj + h.c.) - /i XI (3.1) 

in the two sublattices (even A and odd B) scheme, defining hi = ai and h^j^i = hi 
for z G v4, as this will be useful in the later treatment in presence of holon pairing that 
distinguishes the two Neel sublattices. 

The two fields defined within the magnetic Brillouin zone (MBZ) and diagonalizing 
Hq are 

^±{k) = a^^±e^H^, (3.2) 

— * 

where k E MBZ and the Fourier transforms and 6^ are periodic outside the MBZ. 
The energy eigenvalues are e±{k) = —fj, + ± where = — 2t[cos(A;a;) + cos{ky)] and 
t'^ = —4:t' cos{kx) cos{ky). 

Both fields have pieces of FS within the MBZ (see Fig. O panel (a) ). ^/'^ has four 
hole-like (FS increases as doping increases) arcs centered at Ki, in the middle of the sides 
of the boundary of the MBZ, while ip^ has electron-like (FS reduces as doping increases) 
Fermi arcs near Qi,i = 1, ...,4, at its vertices. In terms of the fields defined in the MBZ 
the holon field, defined in the whole BZ, turns out to be 
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Mk), if ^eMBZ, 
^ \ tfj_{k-2K^), if k^MBZ,ke BZ. 

with {v = 1,2,3,4) chosen to keep the argument of i/)- inside the MBZ (where it 
is defined). 

Now we exploit the two reciprocal primitive vectors Q± = (±7r, tt) to translate the 3rd 
and 4th quadrants respectively in order to obtain a upper rectangular zone D = {k^ & 
[— 7r,7r], ky G [0,7r]} equivalent to the MBZ, as done in the PG (Fig. [2]panel (b)). Notice 
that each translation exchanges the ± index of the fields, because of the minus sign due 
to the exponential factor in eq. (13. 2p . and of the eigenvalues because %_,_q^ = — tg. 

In the presence of holon pairing the field ip^ extended by periodicity, restricted to the 
right quadrant, Dr, and to the the left quadrant, Dl, in D is denoted by \I'+,a, a = R,L. 
It has a good continuum limit because it has a closed hole-like Fermi surface centered at 
Ki in Df> and at K2 in Dl, as discussed in a slightly different framework in [19]. 

Even if mostly neglected in what follows, the field ip- extended by periodicity, and 
then denoted by \E'_, has a small closed electron-like Fermi surface around the point Qi, 
see Fig. [2] panel (b). It is however convenient to restrict also \I'_ to D^i and Dl to treat it 
consistently with ^P+^a, these restrictions are denoted by \E'_,q,, a = R,L. We now discuss 
the holon attraction arising from (12.61) . This interaction is able to distinguish between 
the two sublattices, therefore, when incoherent holon pairs appear the description in the 
MBZ is appropriate. However the pieces of the original FS of the Nearly Free Electron 
description will deform to reach orthogonally the boundary of the MBZ as a consequence 
of the interaction between the previously defined segments of Fermi surface inside and 
outside the MBZ, and the hole- like and electron-like FS discussed above split-off, the fields 
^ having a closed FS. In the BCS treatment adopted in the following the deformation of 
the FS is concentrated near the boundary of the MBZ, where hole-like and electron-like 
FS are closer. Finally we point out the relevance in our discussion of the negative next- 
nearest-neighbor hopping term t' which bends the FS of holons and allows BCS pairing 
of quasi-particles having a good continuum limit in the SM region. In the simple t-J 
model (where t' = 0), the FS in the SM does not cross the boundaries of the MBZ and 
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our approach to the phase fluctuations of the holon pair order parameter presumably is 
possible only in the PG region where closed Fermi surfaces centered at |(±vr, ivr) arise. 
Therefore we guess that in the pure t-J model the formation of holon pairs only appears 
close to T* and not at the naively deduced Tp/^, which instead can give a rough estimate 
of the crossover in the full t-t'-J model. 

The structure of FS for holons discussed above will be qualitatively inherited by phys- 
ical holes, via gauge coupling to spinous discussed later. Such structure bares some 
resemblance with that appearing in the spin-density wave approaches [19,20], but in our 
case its origin is the holon-holon pairing interaction distinguishing the two sublattices, 
not directly the standard AF interaction, although since the holon pairing originates from 
the J term it is still of AF origin. Coexisting hole- and electron-like FS at intermediate 
dopings appear also in large-U treatments of the 2D Hubbard model in terms of an effec- 
tive t-t'-t"-J model [21]. In both cases going from low to high doping level, one first find 
small pockets like in our PG, then coexistence, like in our region between T* and Tp^ and 
finally a large FS, as above Tph in the SM. 

The attractive interaction between holons in different sublattices, hence with opposite 
vortex chirality, of eq. (12. 6p is treated at large scales in analogy with the treatment in 
the PG for the t-J model. Since not all vortices form pairs, a finite screening effect 
persists and the gas of vortices still have a finite correlation length [22], which we denote 
by ^ ~ {Jkp)~^^'^, where kp is the average Fermi momenta of the corresponding closed 
FS. We keep track of the screening effect by replacing in the long wavelength limit A^^ 
in eq. (12. 6p by an effective potential Kff(0 ~ 9/{q'^ + ^~'^)- Hence we approximate the 
interaction hamiltonian at large scales as: 

Hi=-Yl V{i-3)a*h]hjai ^ ^eff(P " q)ap*_^-,-ag- (3.4) 
hJ p,q 

We perform the translation discussed above, add right R and left L labels {a = R, L = 

1,2) to distinguish the two hole-like FS in D and we measure momenta from Kji = Ki 

and Kl = K2 respectively. We neglect the interaction between R and L sectors, except 

for the effect of deformation near the MBZ boundary discussed above and taken into 
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account only phenomenologically, and we adopt the BCS approximation for the R and L 
sectors defining the order parameter 



(3.5) 



For the fields (^1^+, ^E'-, ^E'l, '^\) the Hamiltonian kernel appears in block-diagonal form: 



I ^ a,k OL,k 




a,k 











A'^ . 

a,k 



-A^ . 

a,k 



\ 



K:k ^'^u-'^s 

/z + 1' + 1^^, J 








(3.6) 



with 



ta,k = h+Q. = -2^N^ - (-1)" Sin Cfc = ^'fc+Q. = -4t'(-l)" Sin A;, sin ky (3.7) 

provided A|^ ^ is odd in k; the symbols ± are omitted here and often in the following, 
when their presence can be obviously understood. The choice of the order parameter with 
lowest angular momentum and hence energy has then p-wave symmetry. We assume that 
this order parameter has + Py-wave symmetry for ^P+^l and — p^^-wave symmetry for 
\E'+,ij, so that 

A;_^^ = A^(|^|)7,(fc), 7,(^) = sinfc,. + (-l)"sinfc, (3.8) 

with A^(|fc|) assumed constant, and denoted A^, near the FS in the BCS approximation. 
From eg. (13. 5 p we get ^ = A^ ^ = A^. Then we obtain for the field h a (i-wave 
symmetry gluing R and L sectors when the full BZ is restored (as we see from its definition, 
the order parameter changes its sign when the 3rd and 4th quadrants are translated). This 
p-wave symmetry in D generates also an s-wave pairing for the electron-like Fermi surface 
of as discussed in a similar situation in [19]. 

The energy spectrum of the quasi-particles described by \1'+,q- has the BCS form E^^ = 
±e , r where 

a.±,k 



i^^U-tk-^'y + \K,k\'■ (3-9) 



a,±,fc \ V a,k a,k 
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Let us briely comment on the PG, where 

= -2t ^J {sin k^y + {sin ky)^. (3.10) 

measuring the momenta from Kr {Kl) in the R (L) sectors of D. In the PG the analog 
of the ^'--field in D has no FS, but there is a strong matrix effect due to the Dirac 
structure of the holon field, eventually leading for the hole to an angle-dependent wave- 
function renormalization constant, Z{kx., ky) fa 1 — (sin kx — (—1)" sin ky) / (-\/2((sin k^Y -\- 
{sinkyYY^"^) in Da, even in absence of holon pairing [3]. For the holon pairing the main 
differences between the SM and the PG are the shape of the FS and the higher value of 
the (average) Fermi momentum of the former w.r.t. the latter {kp ~ S), with an induced 
difference for the modulus of the holon-pair field. The Fermi velocity instead is roughly 
the same. 

4 Integration of high energy modes 

In this section we turn to the path-integral formalism and derive the low energy effective 
action which describes the pairing process in the normal state as we lower the temperature 
until the superconducting transition is reached. 

We start from the spinon sector. The relevant action is the non-linear cr-model action 
with a mass term discussed in section 2. This allows us to safely integrate out spinon 
degrees of freedom obtaining an effective action for the slave particle gauge field A. By 
gauge invariance the leading order is a Maxwell-like action: 



37rms 7[o,/3]xr2 

For the holon sector, in the BCS approximation discussed in section 3, the holon 
is gapless only at the 4 nodal points of A|. However in a large-scale gauge- invariant 
treatment whereas one can keep constant the modulus of the order parameter A'' near 
the FS as in BCS, we must include its spatially dependent phase, which we denote by 
0^(x) . This is done setting: 

A^(f,xo) = AV-^^^^'^o), (4.2) 
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(A precise procedure to go from the lattice to the continuum phase field is discussed in 
[23].) The effects of (j)^{x) on holons is non-trivial, as first suggested in a different setting 
in [24], and will be discussed here in detail. We generalize the BCS interaction between 
the holon pairs of section 3 and the FL quasi-particles by allowing phase fluctuations of 
A^. The holon action for \E'+ near the FS in terms of Matsubara frequencies fco then 
reads: 

^^'^(v^^,, A^) = Y.[ik, - E{k) + ^^Wc.^^^.,k (4.3) 

a,k 

ct,k,p 

where we put k = (k^, k). Linearizing the dispersion we set Eik) — ji ^ f^d^l — kpiO)) 
near the FS, with 9 an angle parametrizing the FS. We notice that in the case of constant 
order parameter A|^(A;) — )• A'^5^(A;), Eq. (14. 3 p reproduces the standard BCS coupling 
discussed in section 3. 

One could then reinsert in the kinetic term the gauge fleld by Peierls subsitution 
and naively one would integrate out the holons, again considering the leading term in 
A. However this is not correct since, as will be shown in section 5, above T^. the phase 
fluctuations have short-range decay, with an energy scale that we denote by and 
therefore only the holon modes with higher energy can be safely integrated out. In 
this way no singular terms due to the integration of gapless excitations arise and both 
dynamics and interactions of the order parameter will be approximately local in space. 
A similar approach was proposed in [25], but directly for the hole and using different 
additional approximations. We assume that the high-frequency integrated modes are 
superconducting modes and this assumption will be proved to be self consistent later. 
The result of integration of these high-energy modes in the presence of the gauge fleld A 
can be deduced from the Anderson- Higgs mechanism: the gap equation (13. 5p for each a 
has a degenerate manifold of solutions for arbitrary phase cf)^, therefore the energy should 
depend only on gradients of 0^. Then by gauge invariance of the ho Ion-gauge system we 
obtain in the continuum limit (we take vp = 1 and summation over repeated = 0, 1, 2 
indices is understood henceforth): 
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S^M, ^) ~ E j - 2^. + 2vrn,«)' (^) (4-4) 

where are positive functions and the integer- valued vector currents n^a{x^t) allow 
0^ to be self-consistently an angle function taking into account the presence of vortices 
in the phase of the holon-pair field. Summation on n^a is understood in the partition 
function. 

The formalism of integer currents can be made precise with a lattice regularization 
which we use in the following. One should finally add the action describing the coupling 
of the low-energy modes of \E' with A and (j)^ derived from ( 14. 3 p with UV cutoff m^. 

5 Propagator of the holon-pair field 

In this section we argue on the basis of a self-consistent argument that the Euclidean 
correlation function of the phase of the holon-pair field because of its coupling to a gapless 
gauge field for T > T^. has a purely exponential decay: 

Ga(x, xo) = (e^I'^''(->o)-0''(o,o)]^ ^ e— (5.1) 

This behaviour was assumed in [24] in a different setting; here we explain the origin 
of this behaviour in our approach. (We neglect here for simplicity the index a, as it turns 
out that is the same for both values of a = R,L.) The relation flS.ip defines the 
holon-pair phase coherence length ^ ^/m^f,. can be thought as the mass of the 
phase (j)^ of the holon pairs field and as the mean distance between phase vortices. 

The important point of Eq. (15. ip is that besides its exponential decay for large |x| 
it does not show any singularity for small |x| (in particular it is a constant at x = 0) as 
it would occur for a true massive field which exhibits Ornstein-Zernike decay with power 
in 3D. It is the above stressed feature that eventually is responsible for the peaks in 
the spectral weight precursors of SC, yielding an "effective" gap in the antinodal region 
for the hole, a mark of pseudogap phenomenology. 

We proceed as follows: we first calculate the Ga assuming as action the one obtained 
in section 4, integrating the high-energy modes of spinous and holons and we show that 
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it is self-consistent with the assumed superconducting-hke behaviour of the high-energy 
modes of the holons. Then in section 6 we calculate perturbatively the effect of (f)^ with 
such propagator on the low-energy holon modes and we show that the derived formula for 
the holon propagator if extended to high energy is consistent with the superconducting 
behaviour previously assumed. 

To prove Eq. (15. ip in the setting described above we evaluate the correlation Gi^{x) in 
the Coulomb gauge diAi = (i = 1, 2 indicates spatial components and x = (x, xq)) with 
a lattice regularization. A gauge-fixing is necessary in view of Elitzur theorem [26] that 
states that without gauge-fixing all non-gauge-invariant correlators vanish. Due to the 
presence of the integer vector currents n in (14.41) one cannot evaluate it perturbatively. 
Therefore we use a construction due to Dirac [27] to circumvent the problem: We introduce 
the gauge invariant field 



where E^{y) = 5l^{yo){Ef{y),E^{y),0). Here Ef is the classical (lattice) 2D electric 
field generated by a unit charge at x, hence satisfying the 2D Gauss law diE^{y) = S'i{y) 
and (^"^ denotes the (i- dimensional Kronecker delta. Since Eq. (15. 2p reduces to e*''^''*^^^ in 
the Coulomb gauge, we can evaluate the correlation Ga{x) estimating the expectation 



"electric" field E^{y) has no temporal component and it is different from zero only in the 
temporal plane yo = xq (see Fig. [3]). 

Taking into account the periodicity of the phase, i.e. the presence of vortices, we apply 
the Poisson summation formula to the lattice regularization of (14. 4p rewriting 



(5.2) 



value (e*['^^*^^') </'^(o)]^ without gauge fixing, because it is gauge- invariant. Notice that the 




(5.3) 



n 



J 



where J^(x) are integer vector currents and 




X 



(5.4) 



Then we have: 




(5.5) 
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where VAD(f)^ denotes the standard measure for lattice fields. Defining J^y{z) a current 
of charge 2 supported on a path from x to y (i.e. d^J^y{z) = 2(5^ — 6y){z)), we can write 

<P\x)-<p\0) = J2JU^)d,<l>''i^) (5.6) 

z 

and performing the functional integration in 0'*, appearing linearly in the action, we 
obtain the constraint d^J^{z) = 2(5^ — Sq){z) in the numerator of eq. (15.51) and the 
constraint d^J^{z) = in the denominator. The final step is to integrate out the slave 
particle gauge field A^. The result at zero temperature is: 

-E.^ -3^E.„[^MW+2i?^5W~2i^O(^)]A-l{^-,«)[J^(^«)+2i^-(^«)-2i^O{«;)] 



- E. -^-^ J: Mz)A^\z-w)Mw) 



where Aj^(2;) is the inverse of the 3-dimensional lattice Laplacian (at finite T we just 
replace it by its finite temperature version and restrict the time axis to [0, 1/T] with 
periodic boundary conditions on the fields). The sum in the numerator is on currents 
starting from the point z = Q and reaching the point z = x and on closed currents, 
while only closed currents are present in the denominator. Hence the denominator can 
be interpreted as the partition function of a gas of current loops interacting via the 3D 
(lattice) Coulomb potential. In the numerator the open current has endpoints where the 
electric current spread out in fixed time planes, as described by E . The currents J can be 
interpreted as the Euclidean worldlines of charge 2 particles. Currents supported on loops 
correspond to worldlines of virtual particle-antiparticle pairs, the open current corresponds 
to the worldline of a particle created at one end of the line and annihilated at the other 
one. A similar problem has been dealt with in [28] and the result can be summarized 
as follows: for sufficiently small the closed currents in the leading approximation just 
weakly renormalize the coefficient of the Maxwell action for A, so one can approximate 
(15.71) setting the denominator to 1 and retaining in the sum of the numerator only a single 
fluctuating open current J with weighting factor 



(5 
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The factor e ^"t^ produces an exponential decay in \x\. The power decay of the cor- 
relation function is then decided by how strongly the current J can fluctuate. Gaussian 
fluctuations would produce a power of however in our case the Coulomb potential 

in 3 (Euclidean) or less dimensions is confining. Every fluctuation away from a current of 
minimal length can be described as the addition of a closed loop to the minimal current. 
Because of confinement such loop produces an exponential factor decaying with its area 
and thus strongly suppressing every deviation from a current of minimal length. This 
forces the fluctuations of the current to be non-gaussian, lying within a thin tube sur- 
rounding the shortest straight path (see Fig. |3]). An approximate evaluation of Eq. (15. 8p 
is given in the Appendix and it reproduces Eq. ( 15. ip with a m,^ depending on T through 
A^(T) (contained in c^). The role of the gauge field A was crucial in obtaining this result. 
Without it the open current in (15. 7p would have gaussian fluctuations and the correlation 
would exhibit the standard Ornstein-Zernike decay of a free massive field. 

In Fig. HJwe show the typical temperature behaviour of A^{T) and m^{T) obtained 
with the approximations described in Appendix . 

6 Holon self energy and self-consistency 

In this section we extend the discussion of the previous one about holon pairing considering 
the scattering of holon quasi-particles against the fluctuations of the phase of the holon- 
pairs field. This section was inspired by [24], however here we give an explicit analytical 
formula interpolating from the Fermi-liquid and the superconducting behaviour of holons 
which differs from the one proposed there and in particular our treatment of inclusion of 
dissipation is completely different from that adopted there. We exploit the action in Eq. 
(14. 3 p and Eq. (15. ip to evaluate, in the continuum limit, the zero temperature self energy 
S((X', k) of holon quasi-particles. 

What follows holds both in the PG region and in the SM region since the FS shape 
is irrelevant; for concreteness we discuss the SM case. The difference between the two 
regions lies in the value of the involved parameters, essentially kp. The only condition we 
need is parity invariance which implies that the quasi-particle excitation energy verifies 
E{k) = E{-k). 
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We consider the holon right + field \1/+,_r, the left is similar, with ■yR{k) — )■ jLik). 
Also the field can be treated along similar lines with obvious changes. In the limit 
|A^| « ep the self energy reads 



k) = \A^f^Rikf27cm^^Iik), (6.1) 



Hk) = I 7K±.-r--2- . . r. - (6-2) 



(27r)3 g2 + .(g^ _ k^) _ E{q - k) ' 
where the q dependence of 7r has been neglected for |g| << \k\ ~ kp » and we 
used Gi\{qo,q) = {St^ui^) / {q^ + m|)^. Considering quasi-particles on a shell of thickness 
2 A << kp about the FS we can linearize the quasi-particle dispersion and assuming 
(f « k"^ we have E{q — k) = vp{\q — k\ — kp) ~ vp{—q- k + 6k), where 6k = \k\ — kp, 
obtaining: 



1 J kl + 6k'^ + - 

^^^) ^ T- irrrjL (6-3) 

47r iko + 6k 

Inserting Eq. (16. 3p into Eq. (16. ip and replacing 6k by E{k) we obtain for the self 
energy 



S(a;,A;) = |A'^|V(fc) 



lUJ 



+ E{k)^ + ml 



(6.4) 



where ^nik)"^ ~ 2fc|, sin^ and 6^: is the angle between the vector k and the nodal 
direction Qi in the right sector. The structure of the self-energy (16.40 is reminiscent of 
the electron Green functions of one- dimensional models with dynamically generated mass, 
when the charge and spin velocity coincide [50]. 

The quasi-particle Green's function of the holon is given by 



G{lo, k) = [icu - E{k) - T.{u, 



1 + |A'^|27^(A;)^ 



1 



/w2+£;(fc)2+m2 



(6.5) 



[iLO - E{k)] 
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This quite complicated expression is our key result. We notice that at zero frequency, 
cu = 0, the only pole of the Green's function G{uj, k) is located at Eiji) = 0, that is at 
the FS. This means that there is always a FS except exactly at = 0, which occurs 
only in the superconducting phase. For the existence of the FS the second term in (16. 4p 
is crucial, without it a gap would appear. The absence of the gap is consistent with the 
fact that in spite of a non- vanishing |A^| no phase transition has occurred. 

To understand the meaning of G{co, k) we analyze the two limits of low {u « m^) 
and high [uo » m^) frequency, which coincide with those in [24]. 

For u « niSi we expand the self energy in powers of — up to the second order and 
the Green's function becomes 

Giuj^k)^-^^^ Z^{k)= J (6.6) 

Thus for low frequencies the effect of holon pairing appears through the wave func- 
tion renormalization weight Z/^{k). The system behaves like a FL with unchanged Fermi 
velocity vp and FS of the Nearly Free Electron approximation, but with a strongly direc- 
tion dependent weight Z/^{k) that heavily suppresses (if |A^| >> m^) quasi-particles in 
antinodal directions reducing the effective FS. For u » m^, we expand the self energy 
in powers ^ and we get the Green's function 

G(lo, k) ~ — — (6.7) 

iu^ + E{ky + lA^^l^-fnik)^ 

which is the quasi-particle Green's function of a p-wave superconductor. Thus for 
high frequencies the holon system behaves like a d-wave superconductor in the MBZ, 
the d-wave being obtained gluing the — Py and p^ + Py behaviours in Dji and D^, 
respectively. 

As soon as > 0, the FS and the FL holon quasi-particle peak of weight Z/^{k) 
appear for low frequencies and the system behaves like a metal. When ^ l^^l; quasi- 
particles begin to scatter strongly with the quanta of the phase of the order parameter 
changing significantly both quasi-particle dispersion and coherence. FL quasi-particles 
modes are suppressed in favor of (i-wave SC modes that become much more coherent 

20 



and relevant. This means that decreasing the value of one drives the behaviour of 
the holon system from metallic towards superconducting. Therefore the energy scale m<^ 
separates FL from SC modes, consistently matching with the cutoff in section 4. In fact, 
fermionic modes of energy higher than were previously assumed to be SC-like and 
had been integrated out to give dynamics to the holon-pair field, see eq. (14. 4p . In this 
way one can justify the assumption on high energy modes made in sections 4 and 5 to 
compute the propagator of the phase field, because the obtained self-energy is consistent 
with the assumption made. 

Performing the analytic continuation iu ^ u + ir] in Eq. (16.51) we obtain the retarded 
Green's function 

G^iu, k) = u + ^, + ECk) 



(w + zr/)2 - E{ky - I A'^|27r(A;)2 



1 - 



and its imaginary part is proportional to the spectral weight A{uj, k) = —^QG^{u, k) of 
the holon. 

Since there are FL modes, the scattering rates for holons, ?7, is dominated by Reizer 
singularity and can be computed as in [29]. Combining (16. 5p with loffe-Larkin rule one 
can show that the reduction of the spectral wight reducing temperature yields a deviation 
from below of the T-linear behaviour of in-plane resistivity, typical of the SM , as discussed 
in [30]. 

Finally we should remember that the continuum fields whose spectral weight has been 
discussed above should eventually be converted again to the original holon field using 
(13.31) . obtaining in D in this approximation 



(r/i)(u;,fc) = {{^\^^^+^n){u,k) + {^l^^^+,L){uj,k)) xikeUBZ) 

+ (^l*_)(w,A?)x(A?^ MBZ), (6.9) 

where x denotes the characteristic function. As a consequence of (16.91) the spectral weight 
of the holon coming from the fields does not have peaks outside the MBZ, i.e. it 
only exhibits Fermi arcs along the original FS, except near the MBZ boundary where 
the FS of the NFE approximation is distorted. The same is true inside the MBZ for 
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the field, whose spectral weight therefore exhibits peaks only outside the MBZ. This 
suppression of the spectral weight in the outer "shadow" FS is consistent due to the 
destruction, produced by the fluctuations of the phase field 4>^, of the long-range order 
appearing in the BCS approximation for the holon pairing. A similar suppression in 
the spin-density- wave approach was found in [31], but there was due to fluctuating local 
anti-ferromagnetic order. In this respect physically more similar to our approach is the 
phenomenological model of [32], where the suppression of the outer part of the FS is also 
due to pairing. However, there mathematically this is realized through the presence of a 
zero in the electron Green function, as in the cluster DMFT treatment of the 2D Hubbard 
model [33], which has no direct counterpart in our holon Green function. As a result, in 
those approaches the spectral gap in the antinodal region is asymmetric w.r.t. the Fermi 
energy in contrast with the symmetric behaviour found in our approach below Tph in the 
SM. 

Assuming the Green's function in eq. (16.81) we plot in Fig. [5] panel (a) the holon 
spectral weight at different points of the FS of in the SM with fixed small t]. The 
spectral weight is symmetric in frequencies, and exhibits three maxima, the FL-like peak 
at u; = and two symmetric SC-like peaks roughly at 



iOscp ~ ±^\A^^\^iR{k)^ + ml- (6.10) 

Approximately (al least for small 77) as long as ^ < j^^rp^^5 the FL peak is higher 
than the SC peaks which become negligible for Uscp » l^^llRik), while approximately 
for — > --T^r-^-^, the SC peaks are higher than the FL one. The size of the area under 
the peaks determines the main behaviour of the system. 

Even if small, the FL peak is always present unless strictly vanishes and in this 
case the spectral weight at the FS does not vanish only because of rj. A key point is the 
direction dependence of the condition on the maxima, due to 7R(fc), which determines 
the effectiveness of the FS: 



vlA'^lMk) < [>]ml, (6.11) 
corresponds to a mainly FL[SC] behaviour, respectively. Indeed near the nodal region the 
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behaviour is always FL-like with an effective FS while near the antinodal directions most 
of the spectral weight is concentrated on the SC-like modes and even if there is a FS, it 
has negligible effects. 

The density of states (DOS) of holons DOShioj) can be obtained by numerical integra- 
tion of the spectral weight A{uj^ k) in momentum space. Fig. [5] panel (b) shows the result 
of this integration in the previously defined momentum shell around the FS for different 
values of m<^, using the imaginary part of the retarded Green's function in Eq. (16.81) with 
fixed Tj. 

We see that decreasing the value of the fiat FL density of states gradually reduces 
for small frequencies and develops a peak at a; ~ Wscp- It is a precursor of the ho Ion SC 
peak. Higher frequencies are not affected and preserve a flat DOS. As expected, for small 

the DOS resembles that of a (i-wave superconductor with well defined SC-like peaks, 
but with a finite value at a; = due to the FL quasi-particles which, until vanishes, 
preserve a FS, even if not very effective. 

As shown in section 7, the behaviour strongly dependent on the value of of the 
holon spectral weight leads to a smoother one for the hole spectral weight, where the 
hole is obtained coupling the holon to the spinon through gauge fluctuations, because this 
coupling introduces an additional scattering rate that broadens the peaks appearing in 
the spectral weight of the holon. 

A brief comment for the PG region: the holon spectral weight is somewhat similar to 
that in the SM, but, as shown in Fig. |5] panel (c), due to the linear dispersion the holon 
DOS in the PG exhibits a background linearly increasing with cj, besides a dip and two 
peaks analogous to those in the SM. 

We conclude noticing that if we had neglected the Maxwell-like slave-particle gauge 
field in the computation of the correlation function of the holon-pair field, Ga in section 
4, the fluctuations of the currents defined in section 5 would have been gaussian and 
we would have obtained the standard massive Ornstein-Zernike behaviour instead of eq. 
(15. ip . In this case, using the result in eq. (l6.3p . we can again evaluate the self-energy and 
the Green function of the quasi-particles. The resulting density of states is still reduced 
decreasing at low frequencies, however the two superconducting peaks observed above 
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Tc are now nearly absent. As we shall see in section 7, the SC-like peaks for the holons 
produce analogous peaks in the density of states of physical holes. Their strong reduction 
appearing for standard fermions, due to the absence of the gauge field, disagrees with 
experimental data of tunneling in the cuprates. This argument points out, although 
indirectly, the physical relevance of the slave-particle gauge field. 

7 Reconstructing the hole 

The Green function G for the hole is given in coordinate space by the product of the 
holon and the spinon propagators, averaged over gauge fluctuations. If we now rein- 
sert the gauge fluctuations by Peierls substitution, since the holon modes with frequency 
\uj\ < basically behave as in a FL, their integration provides a Reizer singularity in the 
gauge propagator, thought with a reduced coefficient in the effective action w.r.t. the one 
appearing in absence of holon pairs. The higher frequency holon modes and the spinous 
yields only a sublcading Maxwell-like correction. To obtain a bound state, however, we 
need to take into account the gauge effect non-perturbatevely. To achieve this goal, of 
course in approximate form, in [34] we applied a kind of eikonal resummation of (trans- 
verse) gauge fluctuations; details can be found in the above reference, here we just outline 
the key ideas involved, to explain later on the modiflcations needed for our case. The re- 
summation is obtained by treating flrst A/^ as an external fleld, expanding the correlation 
function in terms of flrst-quantization Feynman paths, then integrating out to obtain 
an interaction between paths which is then treated in the eikonal approximation. Finally 
a Fourier transform is performed to get the retarded correlation function. Whereas the 
path-representation is straightforward for spinous since they are massive, implementing 
it on the holons that have a FS causes problems, due to the presence of an arbitrary 
number of closed fermion wordhnes, describing the contributions of holons in the flnite- 
density ground state. To overcome this difficulty, we apply a dimensional reduction by 
means of the tomographic decomposition introduced by Luther and Haldane [35] to the 
^ fields and then gauge it by minimal coupling. To treat the low-energy holon degrees 
of freedom we choose a slice of thickness A = kp/X, X » 1 in momentum space around 
the FS of the holon. We decompose the slice in approximately square sectors; each sector 
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corresponds to a quasi-particle field in the sense of Gallavotti- Sliankar renormalization 
[36], [37]. Eacli sector is characterized by a unit vector n{6), pointing from the center of 
the FS to the center of the box, labelled by the angle 6 between this direction and the kx 
axis. The contribution of each sector can be viewed approximately as arising from a quasi 
ID chiral fermion; this avoids the finite FS problem for the path- integral representation 
discussed above. The paths appearing for a sector are straight lines directed along the 
Fermi momenta of the sector, with small Gaussian transverse fluctuations. We apply the 
eikonal resummation of gauge fluctuations to the composition of the paths for spinous 
and those arising from holons in each sector. Fortunately the contribution coming from 
the gaussian fluctuations in a sector turn out to be subleading, so in a sector one can 
use Gorkov approximation for holons, i.e. the product of the free Green function time a 
straight phase factor exp(i JJ A^j^dx^). Let us turn to our case. To derive a gauged version 
of Luther-Haldane decomposition for our holon Green function (16. 5 p appears too compli- 
cate, so we proceeds as follows: we approximate (16. 5 p with its limiting behaviours: (16. 6p 
for u < m^, (16. 7p for u > m^. In both cases one can perform the Luther-Haldane decom- 
position, in the second case one should combine a sector with its opposite to construct a 
massive quasi ID fermion [37] and we evaluate the large scale behaviour in approximate 
form. For each sector the Gorkov approximation gives again the leading term. The final 
result of the resummation of the gauge interaction can then be approximately described 
for momenta near the FS by four effects modifying the holon Green function (16. 5p : 

1) the most important one is that the hole inherits from the scattering of gauge 
fluctuations against the spinon, exactly as discussed above in the PG and the SM, a 
renormalized scattering rate, given for T/(X'>>lor<<lby 

r(r,w) ^ {Max{T,uj)KF/t)Q{T,io)/ml (7.1) 

in the SM, where 

Q(T, u) ^ (Max(T, (7.2) 

is the Reizer momentum, and 

r(T, oo) ^ (Max(T, u)KF/t)/ms (7.3) 
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in the PG. In both cases we use the numerical parameters adopted also in refs. [30] and 
[34], respectively. The broadening due to T{T,u) strongly suppresses the contribution of 
the FL holon peak away from the diagonal of the BZ. Here Kp is the (average) Fermi 
momenta of the FS in absence of holon pairing, because the contributions to the gauge 
action come from all the pieces of the FS of the holon h. In the computation we use 
an interpolation between the two limiting behaviours in (17. 2p and fl7.ip . fl7.3p to avoid 
non-smooth behaviour in the plots. 

2) a wave-function renormalization constant Z{T,uj) = {KpmsQ{T,uj)Y^'^ 

3) a renormalization of the chemical potential of the holon by adding to it; we 
denote by E^,{k) the (renormalized) hole dispersion 

4) a sum of two copies of the Green functions obtained as above with Em{k) shifted 
by 2ms. 

The effect 3) is due to the fact that the spinon mass mg naively would produce a gap 
for the hole. As discussed in the PG [34], since the gauge symmetry is unbroken, the 
physical hole must have a Fermi surface and this is achieved by a renormalization of the 
(bare) Fermi momentum which cancel the term m^, guaranteeing a pole at w = in the 
hole Green function. 

The effect 4) is due to the absolute value of in the "relativistic"-like propagator of 
AF spinon, or equivalently to its double-branch dispersion. This phenomenon is absent in 
the standard slave-boson theory because the bosonic holon has " non-relativistic" single- 
branch dispersion. This produces a second term for the hole Green function. As a 
consequence of the shift in the chemical potential discussed above, however the pole in 
the second Green function is shifted with respect to the Fermi energy by 2ms, which is 
approximately the mass of the magnon resonance [3]. One can view the second term in 
the Green function as describing an electron- magnon resonance, due to the "relativistic" 
nature of the spinon. A phenomenologically similar effect appears in the spin-fermion 
approach [38]. As explained in a slightly different framework in [39], this kind of particle- 
hole symmetry breaking is consistent with the original Gutzwiller projection of the t-J 
model. However in our case this asymmetry comes from the coherent contribution and it 
adds to the one arising from the incoherent contribution not considered above. Finally 
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we derive the retarded correlation for a \1/ field by taking the complex conjugation for u 
negative. The resulting retarded correlation near the renormalized FS at an angle 6 is 
given by 

u + ir{T,u) + E<!,{k) 



u + iT{T, a;))2 - E<i,{ky - Al{T) sm{ey S{T, co) 
u + E^{k) + 2m, + ir(T, uo) 



+ 



where 



{uj + iT{T,uo)f - {E^{k) + 2m,)2 - A2(r) sm{dYS'{T,uj) 



SiT,u) = l- 



(7.4) 



(m^(T)2 + E^{ky -{uj + tT{T, a;))2)i/2 
S'{T, 00) = 1 ^ "^^^ . (7.5) 

{m^{Ty + {E^{k) + 2m,y-{u + ir{T,u)yy/^ 

Features of the resulting spectral weight are shown in Fig. [6] panels (a) and (b), where 
the asymmetry and the FL/SC-like crossover with angle and temperature inherited from 
the holon Green function are evident. 

Finally we have to write the full hole retarded Green function combining together 
using (16. 9p those originated from the various \E'-fields considered above, obtaining a Green 
function that exhibits a FS close to the one appearing in the SM without holon pairing, but 
with strongly modified spectral weight. This Green functions qualitatively inherits from 
the holon Green function the pseudogap features discussed in section 6, but with a much 
stronger scattering rate inherited from the spinon scattering against gauge fluctuations. 
These pseudogap features are now directly observable in experiments as discussed in 
section 8, and, contrary to the standard approaches based upon preformed pairs, they 
appear even if there isn't a gas of preformed hole pairs, since the spin degrees of freedom 
are not paired. The term S due to 7^ is crucial to get a FS for the hole even in 
the formal limit of vanishing scattering rate F. If one sets rrifp = 0, hence S = 1, (and 
neglect Z and the second term) in ( \7A\\ the structure of the Green function is of the kind 
discussed in [8]. The F of eqs. (17. 1117. 3p provides the T-dependent pair-breaking term 
introduced phenomenologically there, so many aspects of the pseudogap phenomenology 
are in common in the two approaches. However the non-trivial S term due to 7^ is 
crucial to get a FS without the gap appearing in [8]. This is consistent with the fact that 
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in our approach the appearance of non- vanishing \A^\ occurs from the "normal" states in 
the PG and SM through a crossover, not a phase transition, which seems consistent with 
the experimental data. 

Let us briefly comment on the situation in the PG. Here there are two factors of 
suppression for the hole spectral weight. There is the suppression outside the MBZ, 
due to angle- dependent wave-function renormalization, Z{6), consequence of the Dirac 
structure of holons, and the suppression of the spectral weight away from the diagonals 
due to holon pairing, which is present in the whole PG region, since T* < Tp^. The 
combination of the two effects leads to a spectral weight concentrated, on the FS, near 
the nodes of the SC phase. A quahtative difference with respect to the SM is that in the 
PG there is essentially no "coherent" spectral weight in the antinodal region, since there 
no FS segments are present. 

Finally let us recall that only the "coherent" term due to the hole "resonance" was 
taken into account in the above discussion. For relative momenta between holon and 
spinon larger than Reizer's ~ Q and/or energies larger than ~ vpQ the gauge attraction 
is unable to bind them and one can treat holon and and spinon approximately as non- 
interacting. Their contribution to the hole Green function is therefore obtained as a 
convolution of their free Green functions, producing an incoherent background growing 
for u > (see section 8). A typical hole spectral weight in the PG is presented in Figl6] 
panel (c). 

8 Comparison with experiments 

Let us now summarize the main results of previous sections useful to derive a formula 
for the intensity of ARPES and tunneling. We have shown that in our approach below a 
crossover temperature Tph the holons start to form a gas of holon pairs as a consequence 
of a long-range attraction between spin-vortices centered on holons located on opposite 
Neeel sublattices. The Green function of the phase of the pairing field, minimally coupled 
to the slave-particle gauge field, exhibits a purely exponential decay, with gap (inverse 
correlation length) decreasing with T. In turn the scattering of the fluctuations of the 
phase field against holons produces in the holon Green function, lowering T and moving 
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away from the diagonals of the BZ, a gradual reduction of the spectral weight for u < 
and simultaneously for < uj firstly the formation and then an increase of two peaks of 
intensity of the spectral weight corresponding to SC-like holon excitations, precursors of 
true superconductivity occurring when = 0. Since the holon-pairing distinguishes the 
two Neel sublattices, when it is present the holon fields are naturally periodic in the MBZ, 
the original FS is then distorted, but only near the MBZ boundary for small modulus of 
the holon-pair field. If one consider the MBZ represented by the two upper quadrants 
of the BZ, the holon exhibit two closed hole-like FS centered at the centers of the two 
quadrants, the corresponding continuum field being denoted by and an electron-like 
FS centered at the corners of the region. For > the slave-particle gauge fluctuations 
exhibit a typical scale, a sort of anomalous skin momenta, Q{T,uj) ~ {Max(T,uj)kpY^^ 
and in a range |w| < vpQo around the holon Fermi energy the gauge field couples spinous 
to holons producing a hole resonance. The main contribution to the scattering rate 
T(T,uj) of the hole is due to scattering of the spinon against gauge fluctuations, with 
r(T,w) ^ (Max(r,w))^/3m2 in the SM and ~ Max(T, u;)/ms in the PG. For momenta in 
the MBZ contained in first quadrant of the BZ, the derived explicit form of the retarded 
correlation for the hole resonance can be found in fl7.4p . where 6 is the angle with vertex 
the center of the quadrant and the angle is measured from the direction corresponding 
to the node in the MBZ of the modulus of the holon-pair field. However experimental 
data are usually plotted in terms of the angle, here denoted by a, with vertex the AF 
wave vector and measured from the direction of the node quoted above (see Fig. [7] panel 
(a)). One can get an approximate analytic expression relating 6' to a within the MBZ, 
by considering the holon FS without pairing as a circle centered at the AF wave vector, 
with a radius preserving the area. One obtain thus for an average Fermi momentum 
kp ^ 0.1(1 + 75)1/2 

sin'^~ , ,,,, . (8.1) 

8.1 ARPES 

In our approach the intensity for the electron measured in ARPES experiments due to 
the coherent contribution (of the above discussed hole resonance) for momenta on the FS 
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is proportional to lmGii,[T, —uj,6]n{u)) where n{uj) denotes the Fermi function, and the 
minus sign in is due to the fact that the ARPES experiments deals with electrons, not 
with holes as in our approach. The decrease of the FL peak in the holon spectral weight 
as we move away from the diagonal and lower the temperature yields a similar behaviour 
for the ARPES intensity, but strongly smoothed out by the scattering of spinous by 
gauge fluctuations. The particle-hole asymmetry discussed in section 7 and due to the 
AF structure of spinous smoothed by the thermal broadening appears as a change of the 
decreasing slope of the spectral weight at large positive u. Although the hole spectral 
weight in our approach is not symmetric in u, we can easily take the symmetrization of its 
positive-cj side. Often one has interpreted as a measure of the "spectral gap" on the FS 
half of the distance between the leading peaks of the symmetrized spectral weigh. One 
assumes zero gap for the case of a single leading peak and these "gapless" portions of 
the FS are called "Fermi arcs". In Fig. [7] panel (b) is presented the above defined "gap" 
as function of the angle a, on the portion of the FS relative to \E'+, defined as in (18. ip . 
Taking the above definition of arcs we see that in our scheme their length decreases with 
doping and lowering the temperature, approximately linearly when |A'*| is approximately 
constant and has a T-linear behaviour. However in our approach the appearance of 
the arcs themselves is just an artifact of their definition, the FS is always entirely present, 
although with reduced spectral weight, up to the SC transition, where the FS becomes 
gapped except at "d-wave" nodes. Therefore we propose that the "spectral gap" seen in 
the SM of cuprates is due to thermal broadening of the SC-like holon peaks, masking the 
FL peak. The Fermi arcs then correspond to the region of the FS where the FL peak 
is unmasked. In some respect ours is the opposite of the proposal in [1], where it was 
suggested that the pseudogap phenomenon in the SM should be interpreted as due to 
thermal broadening masking the underlying physical "d-wave" gap. 

In the PG the phenomenology is similar, with three main differences: 1) There is a 
contribution for the FS of the hole corresponding to the holon FS even outside the MBZ 
due to the Dirac structure of holons. However, since Tp^ > T*, when we enter in the PG 
region from the SM the pairing is already active and this implies that the spectral weight 
is already strongly peaked near the node in the MBZ. 2) The interval of energy where is 



30 



effective the attraction between spinon and holon mediated by the slave-particle gauge 
field is much shorter in the PG than in the SM. This region is in fact triggered by the 
Reizer momentum that in the PG is approximately 5"^^^ times smaller than the one in the 
SM. This implies that the contribution of the incoherent part (describing the high energy- 
momentum contribution of the uncoupled spinon-holon) starts at much smaller energy, 
measured from the Fermi energy, thus appearing in the range of energies considered here 
for strong underdoping, in contrast with the SM region. For low T the leading incoherent 
contribution to the hole spectral weight on the FS at energy uj can be approximately 
estimated as DOSh{(^ — ms)/2ms if w > , where DOSh is the density of states of the 
holon. Due to the linear dispersion of the holon in the PG, the DOSh is approximately 
linear in its argument and the incoherent component of the hole gives a contribution to the 
spectral weight at positive u for the hole, growing approximately linearly in uj starting 
from m<( . This particle-hole asymmetric contribution is of the kind discussed in [39]. 
Tentatively we conjecture that the increase of the symmetrized spectral weight observed 
in sufficiently underdoped cuprates moving away from the = peak (see e.g. [40]) 
have (at least partially) this origin. 3) Due to the smallness of kp w.r.t. A'* for strong 
underdoping (except near T*) one finds that is almost 0, hence phenomenologically 
the spectral weight is similar to that of a d-wave superconductor with a large thermal 
broadening due to F, somehow agreeing in the PG with the proposal of [1]. 

Figures [8] show the symmetrized spectral weight in the SM on the FS at 6* = 0.45 for 
different temperatures. The insets present experimental data from [41] exhibiting a good 
qualitative agreement in terms of angle and temperature dependence with the theoretically 
derived behaviour. In particular we see in the data that the dip width is reduced moving 
from underdoped to overdoped samples. We also notice that in underdoped samples, 
described theoretically by the PG, the position of the SC-like peaks is almost constant 
with T, whereas in optimally-overdoped samples, described by the SM, the peaks are 
diverging lowering T, features well reproduced by our derivation. Furthermore, as shown 
in panel (c), we found an approximately linear growth with T, in some range, of the 
spectral weight near the antinode on the FS, due to a saddle point in T. (In the present 
treatment doesn't appear the contribution attributed to pair formation in [41] which 
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would arise from the formation of spinon pairs, not considered in the present work.) The 
dependence on the FS angle a of the symmetrized spectral weight is pointed out in Fig. 
ini panel (a) for SM and (c) for PG and it exhibits a good qualitative agreement with 
experimental data from [41] for SM (panel (b)) and from [40] for PG (panel (d)). 

8.2 Tunneling 

The (SIN) tunneling intensity is usually proportional to the density of states of the elec- 
tron. In the cuprates it might be modified by "matrix element" effects, in particular 
due to orbitals not taken into account by the t-t'-J model and producing a particle-hole 
asymmetric term, but close to u; = such perturbing effects are believed not to produce 
significant qualitative changes [42] and we attempt to compare our theoretical results with 
experimental data. To extract some plot to compare with experiments after approximat- 
ing the holon FS with a circular one with Fermi momentum kp to simplify the calculation, 
we integrate the "coherent" part of spectral weight on momenta in a strip of width ^ Q 
around the FS, to pick up the "coherent" contribution and we have checked numerically 
that almost all the spectral weight fall in this strip. For sufficiently small u even this very 
rough approximation reveals some features of DOS that well compare with the data in 
cuprates. 

We present in Fig. [10] (a) the results for the electron DOS obtained lowering the 
temperature at fixed doping. The DOS exhibits even at high T a broad maximum at 
w = 0, due to the strong non- FL dependence on uj of the scattering rate, furthermore 
below Tph a dip near a; = starts to appear. Both features appear consistent with 
tunneling experiments [43], [44] in slightly overdoped cuprates, that should be described 
by the SM. In a normalized DOS, obtained dividing by the DOS at Tph, (see fig. [10] (c) ) 
the broad maximum disappears. Our DOS is qualitatively similar to a phenomenological 
proposal in [45]. At fixed temperature as the doping increases the depth of the dip 
decreases. 

A final comment on the PG, here the phenomenology is similar except that the in- 
coherent contribution discussed at the end of previous subsection produces an additional 
background growing with energy away from a; = in the electron side. We conjecture that 
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the incoherent contribution is at least partially responsible for the asymmetric behaviour 
seen in SIN tunneling experiments on sufficiently underdoped samples (see e.g. [46]). 



9 Appendix: Computation of Eq. (15.81) 



Following the discussion in section 5 we estimate Eq. (15. 8p in the particular case x = 
{0,Xo), retaining only the straight current, and, to simplify the computation, in the con- 
tinuum limit with ultraviolet cutoff the lattice spacing, e. The open current has only the 
temporal component, Jfj,{z) = 26^o6'^{z)Q{xo — zo), where 9 is the Heaviside step function, 
and, because i?^ has no temporal component, we can decouple from in the term 
with the inverse Laplacian in Eq. (15.71) and factorize two contributions to Ga{x)- The 
first one is due to the straight current J^{z) and the second one is due to the "electric" 
field Ef^: 



where 



27r3 



G'A(a;o,0) ^ e 



(9.1) 



1 / ,2 ,2 ^ ' V 

a xd y 



D{x,) = I £zd'w[^E; - 2El]{z)A,\z - w)[2E; - 2El]{w) = 

'CiXo + C'\n{xo).{9.2) 



Eq. (19.21) follows from the continuum limit definitions of E^{z) = 6{zq — xq) |^^^| {zi 
xi, Z2 — X2, 0) and of the inverse 3-D Laplacian A^^(z) 



with 6 as UV cutoff. 



The last step holds in the limit Xq » 1 » e and Ci,C' positive constants. 

The important point is that the logarithmic contribution in Eq. (19. 2 p is exactly elim- 
inated from the integral in the first exponential of Eq. (19. ip due to the straight current. 



/ dt dt'A^\t-t',6) = dt dt' , ^ C2X0 - C \n{xo) (9.3) 

Jo Jo Jo Jo ^J{t-t'Y + e'^ 

with C2 a positive constant and numerically C3 = Ci + C2 ~ 25. In conclusion Gi^{x) 

turns out to have a purely exponential decay in |x| and from the above discussion, our 

estimate of the inverse coherence length for small A'* is m^f, ~ 2c^ + ^^^rns- 
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Actually with the more sophisticated methods of the excitation expansions of [48], one 
obtains for small A'^ the estimate 

m^^ —-{ - P{c^) (9.4) 

where P is a polynomial taking into account the fluctuations of the current J^{x)] 
in an extrapolation at lower temperatures it ensures the right behaviour when the SC 
temperature for the holon system, Tch, is approached and tends to zero. It should be 
stressed that T^h is not the physical SC temperature T^, because in the physical system 
also the spinous contribute to the SC transition through the formation of RVB pairs, so 
that Tch ^ Tc. 

In BCS approximation Ci = C2 = c is proportional to the superfluid density, ~ (A^)^ 
for small A'^ , plus a correction ~ m| and Cq is approximately proportional to kp/v^, 
therefore for sufficiently small A^ and kp the above computations are justified. However, 
to compute using the above procedure the explicit dependence of on T needed in the 
following would be a formidable task. We expect and we have partially verified numerically 
that the physical results depend only weakly on the explicit form of this dependence and 
in the computations we adopt the following drastic simplifications: First we replace v,/, 
hj vp, this would reproduce the correct order of magnitude for an s-wave pairing and 
according to the results in [49] it is not unreasonable even for a d-wave pairing; anyway 
we have verified that numerically the precise value is essentially irrelevant in the range of 
temperatures considered here. Second we keep only the linear term of P, with a coefficient 
chosen exactly to ensure the vanishing of at the (unphysical) holon superconducting 
temperature Tch , assuming that all the temperature dependence is appearing in A^(T). 
We thus derive the following expression for m^{T) used in the plots: 

m^(T) = {-A'^iTchf + {A'^iTchY + 4A\A\Tchf - A'^{Tf)y/')/{2A). (9.5) 

In (19. 5p A = {2kF)~^ + (757r/2ms) and (away from T = 0) we fit numerically A''(T) 
by the approximate expression A'^(T)/A^(0) ^ [1 - e2-5(i-rph/^))]i/2 ^nd [2] A^(0) ^ 
1.6Z(5)A;y with Z{6) = {I - 26){{1 + m^y/^ - m,). 

Furthermore since the holon system is not in the weak BCS region, Tph » Tch and 
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we just set Tch ~ 0. 
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Figure 1: Schematic behaviour of the crossovers T*, Tp^ and Tpg in the 6-T phase diagram. 
Concentric oriented circles (red onhne) denote the charge flux, oriented filled disks (blue 
online) denote the spin vortices, straight lines their pairing attraction, wavy lines the 
gauge interaction, arrows the spinous. 




Figure 2: Panel (a) represents the FS of the t-t' model (t=0.31eV, |t'|=0.08eV, /i=0.15eV) 
defined in the whole BZ (coinciding with that of the t-t' -J model in the MFA discussed in 
the text in the SM); the holes fill the shaded region. The boundaries of the MBZ are also 
shown. Panel (b) represents the upper rectangular part of the BZ equivalent to the MBZ 
after the translation of the 3rd and 4th quadrants with the original FS (blue lines online) 
and the translated one (red lines online); dashed lines are qualitatively the closed FS of 
\E'+ (centered in Ki and K2) and ^ - (centered in Qi)in the presence of holon-pairing. 
Axis are measured in units of tt. 
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Figure 3: The straight temporal current J^{z) starting from the point z = 0, ending at 
the point z = x and joining the two orthogonal (constant ^o) planes which are the support 
of the electric fields E^{z) and E^{z). Suppressed fluctuations of the straight current are 
also shown. 
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Figure 4: A^{T) and m^(T) for 6 = 0.18 
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Figure 5: Panel (a) shows the ho Ion spectral weight (imaginary part of Eq. (16.81) ) at 
k = kp versus energy and the FS angle 6, in units of vr, in the SM at fixed r] = 0.05J; the 
related holon DOS for different values of T is shown in panel (b). Panel (c) shows the 
ho Ion DOS in the PG for different values of T at the same rj. 
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Figure 6: Panels (a,b) show the hole spectral weight in the SM (imaginary part of Eq. 
f l7.4p )and panel (c) in the PG at /c = fci? as a function of the FS angle 6 (in units of vr) 
and energy. Holes fill the positive energy side. 
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Figure 7: Panel (a) shows the definition of the FS angles 6 and a in the first quadrant of 
the BZ. Panel (b) shows the spectral gap (half the distance between the SC peaks in the 
symmetrized spectral weight) versus FS angle a measured in units of n. In the inset data 
from [47]. 
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Figure 8: Panels (a) and (b) show the symmetrized spectral weight at the "antinodal" FS 
direction 9 = 0.45 in the SM. The curves are subsequently shifted upward by 0.08 with 
decreasing T . In the insets experimental data from [41]. Panel (c) shows the spectral 
weight at the Fermi energy near the antinode for different doping values, renormalized by 
subtracting its value at Tph. In the inset data from [41]. 
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Figure 9: Panel (a,c) shows the symmetrized spectral weight as a function of the FS 
direction a in the SM and in the PG, respectively . Panel (b) shows the experimental 
data for 5^2212 [41] and panel (d) for La2-xSr^CuO^ [40]. 




Figure 10: (a) The electron DOS in the SM . (c) The same electron DOS renormalized 
by subtracting its value at Tyh at different temperatures. The curves are subsequently 
shifted up of 0.04 with decreasing T for clarity. The experimental data (b) and (d) are 
from [44]. 
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